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ABSTRACT 

We interpret the large variety of redshift distributions of galaxies found by far-infrared and (sub-)millimeter deep surveys depending 
on their depth and wavelength using the Bethermin et al. (2012) phenomenological model of galaxy evolution. This model reproduces 
without any new parameter tuning the observed redshift distributions from 100 pm to 1.4 mm, and especially the increase of the median 
redshift with survey wavelength. This median redshift varies also significantly with the depth of the surveys, and deeper surveys do 
necessarily not probe higher redshifts. Paradoxically, at fixed wavelength and flux limit, the lensed sources are not always at higher 
redshift. We found that the higher redshift of 1.4mm-selected south pole telescope (SPT) sources compared to other SMG surveys 
is not only caused by the lensing selection, but also by the longer wavelength. This SPT sample is expected to be dominated by a 
population of lensed main-sequence galaxies and a minor contribution (~10%) of unlensed extreme starbursts. 
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1. Introduction 


The determination of the star formation history in the Universe 
is a key challenge for unde rstanding the evolution of galaxies 
iMadau & DickinsCTl2014l) . Two decades ago, the first deep 
submillimeter surveys (850 pm) reve aled a population o f dusty, 
strongly star-forming galaxies (e.g.. ISmail et alJll997h . which 
were missed by optical surveys. This showed the importance of 
submillimeter observations to draw a complete pi cture of the star 
forma tion activity in the high-redshift Universe. IChanman et al.l 
ll2005l) found that the median redshift of these submillimeter 
galaxies (SMGs) is 2.3. However, the identification of the op¬ 
tical counterparts was non trivial because of the large beam of 
single-dish submillimeter telescopes (~15-20"). Their precise 
position was determined by radio interferometry and their red¬ 
shift measured using optical spectroscopy. This method was thus 
potentially biased against the highest redshift objects and galax¬ 
ies in the 1.4<z<2 redshift desert. 


A decade later, I Vieira et"^ (l201(l 1201 3l) identified a popu¬ 
lation of strongly-lensed dusty star-forming galaxies (DSFGs) 
in the 1.4 mm South Pole Telescope (SPT) survey. Using 
ALMA to derive directly spectroscopic re dshifts by ta r geting 
CO-transitions at millimeter wavelengthSjJWeiB^e^aH (2013) 
found that their median redshift is 3.5. ISimpson et al.l ( 2014 ) 
used millimeter interferometry with a subarcsec resolution to di¬ 
rectly confirm the optical counterparts of SMGs and measure d a 
median photometric redshift similar to lChapman et alJ (l2005h . Is 
this difference of median redshift caused by the different wave¬ 
lengths or the lensing e ffect? In this letter, we d iscuss these se¬ 
lection effects using the iBethermin et al.l (1201 2al. hereafter B12) 
model, that reproduces the redshift distribution of sources de¬ 


tected at all wavelengths from 100^(m to 2 mm (Sect.|2]). Us¬ 
ing a simplified version of this model, we explain why sources 
selected at longer wavelengths have a higher median redshift 
(Sect.O. We then show how the redshift distribution of DSFGs 
is affected by the flux and lensing selection (Sect. m. In Sect.|5] 
we discuss in more detail the SPT-selected population of galax¬ 
ies in the light of our findings of the two preceding sections. We 
define Lir as the bolometric energy emitted by a galaxy between 
8 and lOOOjum. 


2. Model and comparison with observations 


The B12 is model is based on the observed evolution of the stel¬ 
lar mass (M*) function and the main-sequence of star-forming 
galaxies, i.e., a tight SFR-M* correlation. At fixed mass, 
the star formation rate (SFR) increases rapidly but smoothly 
with redshift. This model also includes a population of star- 
bursts with a strong excess of sSFR (^SFR /M.^.) that contribute s 
15% of the star formation density at z >l dSarg ent et al .]l2012l) . 
We used two different SEDs based on iMagdis et al. ( 20121) for 
main-sequence and starburst galaxies. The dust of the main- 
sequence templates is warmer at higher redshifts. The contribu- 
tion of strongly-lensed ga laxies is computed using the model of 
iHezaveh & Holde3(l201 lb . The source size can affect the lensing 
amplification, but an extreme size evolution by a factor of 12.5 
between z=2.5 and z-6 would be required to m atch the 850jum 


and SPT redshift distributions (WeiB et aT] 


vation s dlkarashi et al.ll201^ISmolcic et al 


20131) . Recent obser- 
2014l:[^mpson et ^ 


12015l) do not suggest size evolution. We thus do not adopt any 
size evolution in our model. In this letter we focus on the 100- 
2000 pm range, since at shorter and larger wavelengths the AGN 
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Fig. 1. Comparison between the prediction of our model (orange 
filled histogram) and the observations (see Table[TJ. The flux cuts indi¬ 
cated to compute the redshift distributions are indicated on the figures. 
The distributions are normalized to have J(dNldz)dz - 1 (expect the 
100yum normalized at 0.5 for clarity). 

have a non-negligible contribution throug h their torus and/or 
synchrotron emissions dProuart et al.ll2014li . 

We test the validity of our model by comparing it with a large 
compilation of redshift distributions listed in Table[T] Overall the 
model reproduces well the observed redshift distributions from 
100yum to 1.4 mm (see Fig.ITJ. At 2 mm (not shown in Fig. [D , 
only 5 of the of the 12 sources found bv IStaguhn et al.l (l2014h 
have known redshifts. Our model predicts a median redshift of 
2.9 assuming of flux cut of 0.24 mJjQ, while they found 2.9+0.4. 
Note that the good agreement between our model and the data 
is not the result of fine tuning the ingredients of our model. We 
can nevertheless observe a small tension at 850yum, where our 
distribution peaks at a redshift too high by Az~0.2. At the flux 
density cuts used for all panels in Fig[T] our model predicts that 
we mainly select main-seq uence galaxies. The smooth evolu¬ 
tion of th e main-sequence (ISchreiber et al.|[2014t) and the mass 
function (filbert et al.ll2013h . as well as the volume effects are 

* The flux density of their faintest source after deboosting. 


Table 1. Summary of data used in Fig.fT] 


Reference 

N 

1 a 

Slim 

mJy 

Method 

Berta et al. (2011) 

5360 

100 

9 

a, b 

Bethermin et al. (2012b) 

2517 

250 

20 

a, b 

Geach et al. (2013) 

60 

450 

5 

b 

Casev et al. (2013) 

78 

450 

13 

b 

Wardlow et ak (2011) 

72 

850 

4 

b 

Simnson et al. (2014) 

77 

850 

4 

b, e 

Chanman et al. (2005) 

73 

850 

3 

c 

Smolcic et al. (2012) 

28 

1100 

1.4 

b, e 

Michalowski et al. (2012) 

95 

1100 

1 

b 

Yun et al. (2012) 

27 

1100 

2 

b 

Weifl et al. (2013) 

23 

1400 

20 

d, f 

Staguhn et al. (2014) 

5 

2000 

0.24 

b, c 


Notes, a) extraction of the sources based on PSF-fitting codes us¬ 
ing short-wavelength priors; h) photometric redshifts; c) optical/near- 
infrared spectroscopic redshift after radio identification; d) millime¬ 
ter spectroscopic redshift; e) identification of the optical/near-infrared 
counterparts using high-resolution (sub-)millimeter data. 


thus suffi cient to explain the r edshift distribution from lOOyum to 
1.4 mm. IZavala et akl (1201 4|) extrapolated successfully the red¬ 
shift distributions from 850 /urn up to 1.4 mm, but did not manage 
to reproduce the 450 fim. They claimed another galaxy popula¬ 
tion is necessary. This could be caused by their assumption of 
a flux-invariant redshift distribution, which is not present in our 
more refined model (see Sect. nil. 


3. Why does selection at longer wavelengths select 
higher-redshift galaxies? 

The shift of the peak of the redshift distribution towards higher 
redshift when wavelength increases (Fig.[T]i can be easily ex¬ 
plained using a simplified version of our model. In this ver¬ 
sion, we neglect the effect of the strong lensing and assume that 
all galaxies at a give n redshift have the a verage main-sequence 
SED provided by the lMagdis et al.1 (1201 2l) template library. The 
left panel of Fig.|2] shows how the Lir limit at which sources 
can be detected evolves with redshift for the flux density cuts 
used in Fig.[T] At z>l, the curve is almost flat at 1.4 mm, but 
increases quickly with increasing redshift at 250yum. This effect 
is caused by the shift of the peak of the dust emission (~100yum 
rest-frame) toward the ob served millimeter wavelengths when 
the redshift increases (e.g., iBlain et al.ll2002l) . 

The redshift distribution can then be deduced from the lumi¬ 
nosity function (LF). The right panel of Fig.|2]shows the cumu¬ 
lative luminosity functions at z=l and z=4 computed from the 
B12 model. These functions are renormalized to directly pro¬ 
vide the number of sources per square degree and redshift inter¬ 
val above a certain Lir cut. The x-axis and y-axis switched with 
respect to the usual representation to better illustrate the connec¬ 
tion with the left panel. At 250 yum, the Lir limit of surveys is 24 
times larger at z=4 than at z=l. Consequently, despite a slightly 
higher density of luminous objects at z=4, we detect 450 times 
fewer objects at z=4, because of the steep slope of the luminosity 
function above the knee. At 1.4 mm, the Lir limits are similar at 
both redshift, but we detect 5 times more objects at z=4, because 
of the evolution of the LL. 
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Bethermin et al.: The influence of selection effects of redshift distribution of SMGs 



redshift dN(>LIR limit)/dz [gal.deg'^] 

Fig. 2. Simplified diagram illustrating the impact of the wavelength on the redshift distribution. The left panel represents the detection limit in 
bol ometric infrar e d lum inosity (Lir) as a function of redshift for a S 250 >20 mjy (purple) and S 1400 > 1 mjy (gold). These limits are computed using 
the lMagdis et al.l i2012h main-sequence SED library. We neglected the scatter and the starburst populations to simplify this diagram. The right 
panel represents the number density of objects per redshift interval (computed from the luminosity function and the geometry of the Universe) as 
a function of the Lir cut (x and y axis are inverted) at z=l (dashed line) and z=4 (dot-dash line). These curves are computed using our model. The 
colored arrows represent the ratio between the number density of z=l and z=4 sources. The black arrow highlights the increase of the luminosity 
and the decrease of the density of the knee of the LF between z= 1 and z=4. 


4. Impact of flux and lensing on the redshift 
distribution 

We also used the standard version of our model (including 
strong lensing, different templates for main-sequence and star- 
burst galaxies, and temperature dispersion) to predict more pre¬ 
cisely how the median redshift of galaxies selected by photomet¬ 
ric surveys changes depending on the wavelength and the flux 
density cut (Fig.O. 

At A >450 /im, we found an abrupt decrease of the median 
redshift around 100mJy (lOmJy at 2mm). This sharp transition 
is caused by the presence of two bumps in the redshift distribu¬ 
tions: local galaxies around the knee of the luminosity function 
and luminous galaxies z~2-3. There are fewer sources at z~0.5-l 
than at z~0.1, because the Lir limit is well above the knee of the 
luminosity function. There are also more sources at z~2-3 than 
at z~0.5-l, because the volume probed is larger and the num- 
ber of luminous galax ies per unit of volume higher (see, e.g., 
iBethermin et al.l 1201 ll) . At high fluxes, the number of nearby 
obje cts decreases following an Eucli dian trend in (see, 

e.g. iPlanck Collaboration et al.ll2013h . At higher redshift, we 
are probing objects above the knee of the luminosity function, 
where the slope is exponential. Consequently, when we observe 
at very high flux, the nearby population always dominates the 
redshift distribution (e.g. the S500 > 100 mJy sou rces), while the 
lensed population has a median red shift of ~2.4 (iNegrello et al.l 
l2010tlGonzalez-Nuevo et al.ll2012l) . 

We now consider the redshift distribution of only the 
strongly-lensed galaxies (dashed line in Fig.|2l. At /1<1.1 mm, 
the lensed galaxies are always at higher redshift than the full 
population. This is expected, since the probability of lensing in¬ 
creases with redshift because of the larger probability to find a 
massi ve galaxy on the line of sight (see, e.g.. lHezaveh & Holdeii 
I 2 OI ll) . However, at /1>1.1 mm, the lensed objects are at a lower 
median redshift in some specific flux range (2-7 mJy at 1.4 mm). 
This can be explained by the lower intrinsic fluxes (<lmJy) 
of the lensed sources. The redshift evolution of the luminosity 
function and Lir limit (Fig.|2]i puts these faint sources at lower 
median redshift than the >1 mJy sources which dominate the un- 
lensed population. In the domain of a few mJy, this effect can 


thus compensate the lensing effect, which biases the redshift dis¬ 
tributions toward higher redshift. 

These results suggests that the best strategy to build sam¬ 
ples of high-redshift galaxies is to perform surveys at the largest 
possible wavelength. However, all objects below z=8 have sim¬ 
ilar 1.4 mm/2 mm colors, because they are all observed in the 
Rayleigh-Jeans regime. The flux at 2 mm is 3 times lower, while 
the sensitivity of instruments is similar for the two wavelengths. 
Consequently, the expected number of detections for the same 
integration time is significantly lower. The higher median red¬ 
shift at 2 mm is mainly caused by the lack of Lir sensitivity at 
z<2 rather than a better efficiency to detect very high redshift 
sources. A compromise has thus to be found between the large 
wavelength and the efficiency of the survey, defined here as the 
number of detections per hour. Furthermore, the risk of contam¬ 
ination by free-free emission increases and source identification 
at longer wavelengths is more difficult due to the larger beam. At 
A > 500 jum, the median redshift decreases, when we go deeper 
than 1 mJy. Deeper surveys do not automatically imply higher 
redshifts. 

5. A closer look at the SPT galaxy population 

The SPT survey found a l arge population of l ensed galaxies us- 
ing a selection at 1.4 mm d Vieira et al.ll2010h . IWeiB et ^ d2013l) 
found that their median redshift is the highest measured to date 
in a simple photometric selection. We will discuss their char¬ 
acteristics in this section. Fig. 0] shows the redshift distribution 
predicted by the B12 model of sources selecte d at 1.4 mm with 
the 20 mJy density flux cut of lWeiB et alJ d2013l) . but also a lower 
(1 mJy) and a higher (lOOmJy) one to illustrate how this influ¬ 
ences the nature of the selected sources. The contribution of 
z< 0.1 is significant for 20mJy and 100 mJy (32% and 70%, re¬ 
spectively), but negligible for 2 mJy (0.9%). Since we are more 
interested in the high-z population, we normalized the redshift 
distribution in Fig.|4]ignoring the objects at z<0.L 

The nature of the high-redshift sources also depends strongly 
on the flux cut. For the 20mJy cut, 87% of sources are lensed 
and 90% of them are classified by the model as "main-sequence". 
Using a mock catalog based on our model, we estimated that 
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Fig. 3. Median redshift of dusty galaxies as a function of the flux cut 
at various wavelengths (see color coding in the plot). The solid lines 
correspond to full samples and dashed lines to strongly-lensed samples. 
The filled dot indicates the limit of one source per In sr (about half of 
the sky is sufficiently clean for extragalactic surveys). The observations 
listed in Table[T]are symbolized by diamonds (full samples) and a square 
(SPT lensed sample). 

their average sSFR is 1.7 times higher than the center of the 
main-sequence, while in our model, we define starbursts as ly¬ 
ing > 4x above. Because of their higher SFR at fixed stellar 
mass, the galaxies slightly above (~lcr) the center of the main 
sequence are easier to detect than objects slightly below. The 
unlensed sources are extreme HyLIRG (Lir >10'^ Lq) starbursts 
with an excess compared to the main sequence larger than a fac¬ 
tor of 4. For lOOmJy, all the high-redshift sources are lensed. 
This is expected, because the minimal intrinsic luminosity nec¬ 
essary to detect a source without lensing, assuming our coldest 
template, is 10*^ Lq, which is unphysical. For a cut of 2 mJy cor¬ 
respo nding to the expected sensitivity of SPT-3g (iBenson et al.l 
12014 . we typically expect to detect unlensed ULIRGs ((Lir) = 
7x10*^ Lq) at z~3.2. At this redshift, this corresponds to the 
knee of the luminosity function. The slope of the counts is thus 
shallower and the contribution of lensed sources is small (2.9%). 

Our model predicts that SPT sources lie only slightly above 
and still well within the scatter of the main sequence. The 
high SFR of these objects is expected to be caused by large gas 
reservoirs rather than a merger-induced, boosted star -formation 
efficiency (ISargent et al.ll2014lBethermin et al.ll2015h . 


6. Conclusion 

The B12 model allows us to understand how the observed red¬ 
shift distributions of DSFGs depend on how they are selected. 
Our main findings are: 

- The B12 model successfully reproduces the redshift distribu¬ 
tions from 100 pm to 1.4 mm without any additional param¬ 
eter tuning. 

- When we select sources at longer wavelength, the median 
redshift of the sources increases. This effect can be easily 
explained considering the Lir limit versus redshift of the sur¬ 
veys and the evolution of the infrared LF. 

- The median redshift of the DSFGs also varies with survey 
depth. The deeper surveys in the (sub-)millimeter surpris¬ 
ingly probe lower redshifts. At A >1.4 mm in specific flux 



Fig. 4. Upper panel: redshift distribution from the B12 model of 
dusty galaxies selected at 1.4 mm for various flux cuts (see color coding 
in the plot). The unlensed galaxies are represented by solid lines and 
the contribution of lensed populations by a dashed line. For each flux 
cut, we normalized the sum of the lensed and unlensed distributions 
considering only z>0.1 sources to allow an easier comparison of high 
redshift distributions. Lower panel: fraction of starbursts. The curves 
are not plotted, where dN/dz<l sr“'. 


intervals, the lensed objects can also be at higher redshifts 
than unlensed sources. 

- The DSFGs selected by SPT are mainly strongly magnified, 
main-sequence galaxies, but 10% of these sources are pre¬ 
dicted to be unlensed HyLlRGs. 
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